Metabolite signatures associated with microRNA miR-143-3p serve as drivers of poor lung function trajectories in childhood asthma.

BACKGROUND: Lung function trajectories (LFTs) have been shown to be an important measure of long-term health in asthma. While there is a growing body of metabolomic studies on asthma status and other phenotypes, there are no prospective studies of the relationship between metabolomics and LFTs or their genomic determinants. METHODS: We utilized ordinal logistic regression to identify plasma metabolite principal components associated with four previously-published LFTs in children from the Childhood Asthma Management Program (CAMP) (n = 660). The top significant metabolite principal component (PCLF) was evaluated in an independent cross-sectional child cohort, the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS) (n = 1151) and evaluated for association with spirometric measures. Using meta-analysis of CAMP and GACRS, we identified associations between PCLF and microRNA, and SNPs in their target genes. Statistical significance was determined using an false discovery rate-adjusted Q-value. FINDINGS: The top metabolite principal component, PCLF, was significantly associated with better LFTs after multiple-testing correction (Q-value = 0.03). PCLF is composed of the urea cycle, caffeine, corticosteroid, carnitine, and potential microbial (secondary bile acid, tryptophan, linoleate, histidine metabolism) metabolites. Higher levels of PCLF were also associated with increases in lung function measures and decreased circulating neutrophil percentage in both CAMP and GACRS. PCLF was also significantly associated with microRNA miR-143-3p, and SNPs in three miR-143-3p target genes; CCZ1 (P-value = 2.6 × 10-5), SLC8A1 (P-value = 3.9 × 10-5); and TENM4 (P-value = 4.9 × 10-5). INTERPRETATION: This study reveals associations between metabolites, miR-143-3p and LFTs in children with asthma, offering insights into asthma physiology and possible interventions to enhance lung function and long-term health. FUNDING: Molecular data for CAMP and GACRS via the Trans-Omics in Precision Medicine (TOPMed) program was supported by the National Heart, Lung, and Blood Institute (NHLBI).

Cluster Analysis of World Trade Center Related Lower Airway Diseases.

ABSTRACT: Introduction: Cluster analysis can classify without a priori assumptions the heterogeneous chronic lower airway diseases found in former workers at the World Trade Center (WTC) disaster site. Methods: We selected the first available chest computed tomography scan with quantitative computed tomography measurements on 311 former WTC workers with complete clinical, and spirometric data from their closest surveillance visit. We performed a nonhierarchical iterative algorithm K-prototype cluster analysis, using gap measure. Results: A five-cluster solution was most satisfactory. Cluster 5 had the healthiest individuals. In cluster 4, smoking was most prevalent and intense but there was scant evidence of respiratory disease. Cluster 3 had symptomatic subjects with reduced forced vital capacity impairment (low FVC). Clusters 1 and 2 had less dyspneic subjects, but more functional and quantitative computed tomography evidence of chronic obstructive pulmonary disease (COPD) in cluster 1, or low FVC in cluster 2. Clusters 1 and 4 had the highest proportion of rapid first-second forced expiratory volume decliners. Conclusions: Cluster analysis confirms low FVC and COPD/pre-COPD as distinctive chronic lower airway disease phenotypes on long-term surveillance of the WTC workers.

Socioeconomic status, diet, and recurrent severe asthma exacerbations in Puerto Rican youth.

Background: Why Puerto Rican youths have higher rates of severe asthma exacerbations (SAEs) than their non-Hispanic White peers is unclear. Objective: We aimed to identify risk factors associated with recurrent SAEs in Puerto Rican youths with asthma. Methods: We performed cross-sectional and longitudinal analyses of recurrent SAEs in 209 Puerto Rican youths with asthma who participated in 2 cross-sectional studies approximately 5.2 years apart: the Puerto Rico Genetics of Asthma and Lifestyle study (visit 1, participants aged 6-14 years) and the Epigenetic Variation and Childhood Asthma in Puerto Ricans study (visit 2, participants aged 9-20 years). Recurrent SAEs were defined as at least 2 SAEs in the previous year. Results: Of the youths in our study, there were 80 (38.3%) and 47 (22.4%) with recurrent SAEs at visit 1 and visit 2, respectively, and 31 participants (14.8%) had persistent recurrent SAEs (ie, recurrent SAEs at both visits). In multivariable analyses, low household income was significantly associated with 2.4 to 12.3 times increased odds of recurrent SAEs in all analyses, with stronger longitudinal associations. Low parental education level, nonprivate or employer-based health insurance, overweight or obesity, residential proximity to a major road, and low or moderate level of outdoor activity were each significantly associated with recurrent SAEs in at least 1 analysis. Further, persistence of low parental numeracy level, low household income, and an unhealthy diet were each associated with persistent recurrent SAEs. Conclusion: In this study of Puerto Rican youths with asthma, persistence of low parental numeracy level, a low household income, and an unhealthy diet were associated with persistent recurrent SAEs. Our findings support policies promoting equity and healthy lifestyles for Puerto Rican children and their families.

Asthma morbidity measures across Black ethnic subgroups.

BACKGROUND: Black adults are disproportionately affected by asthma and are often considered a homogeneous group in research studies despite cultural and ancestral differences. OBJECTIVE: We sought to determine if asthma morbidity differs across adults in Black ethnic subgroups. METHODS: Adults with moderate-severe asthma were recruited across the continental United States and Puerto Rico for the PREPARE (PeRson EmPowered Asthma RElief) trial. Using self-identifications, we categorized multiethnic Black (ME/B) participants (n = 226) as Black Latinx participants (n = 146) or Caribbean, continental African, or other Black participants (n = 80). African American (AA/B) participants (n = 518) were categorized as Black participants who identified their ethnicity as being American. Baseline characteristics and retrospective asthma morbidity measures (self-reported exacerbations requiring systemic corticosteroids [SCs], emergency department/urgent care [ED/UC] visits, hospitalizations) were compared across subgroups using multivariable regression. RESULTS: Compared with AA/B participants, ME/B participants were more likely to be younger, residing in the US Northeast, and Spanish speaking and to have lower body mass index, health literacy, and <1 comorbidity, but higher blood eosinophil counts. In a multivariable analysis, ME/B participants were significantly more likely to have ED/UC visits (incidence rate ratio [IRR] = 1.34, 95% CI = 1.04-1.72) and SC use (IRR = 1.27, 95% CI = 1.00-1.62) for asthma than AA/B participants. Of the ME/B subgroups, Puerto Rican Black Latinx participants (n = 120) were significantly more likely to have ED/UC visits (IRR = 1.64, 95% CI = 1.22-2.21) and SC use for asthma (IRR = 1.43, 95% CI = 1.06-1.92) than AA/B participants. There were no significant differences in hospitalizations for asthma among subgroups. CONCLUSIONS: ME/B adults, specifically Puerto Rican Black Latinx adults, have higher risk of ED/UC visits and SC use for asthma than other Black subgroups.

Nasal epithelial gene expression and total IgE in children and adolescents with asthma.

BACKGROUND: Little is known about nasal epithelial gene expression and total IgE in youth. OBJECTIVE: We aimed to identify genes whose nasal epithelial expression differs by total IgE in youth, and group them into modules that could be mapped to airway epithelial cell types. METHODS: We conducted a transcriptome-wide association study of total IgE in 469 Puerto Ricans aged 9 to 20 years who participated in the Epigenetic Variation and Childhood Asthma in Puerto Ricans study, separately in all subjects and in those with asthma. We then attempted to replicate top findings for each analysis using data from 3 cohorts. Genes with a Benjamini-Hochberg-adjusted P value of less than .05 in the Epigenetic Variation and Childhood Asthma in Puerto Ricans study and a P value of less than .05 in the same direction of association in 1 or more replication cohort were considered differentially expressed genes (DEGs). DEGs for total IgE in subjects with asthma were further dissected into gene modules using coexpression analysis, and such modules were mapped to specific cell types in airway epithelia using public single-cell RNA-sequencing data. RESULTS: A higher number of DEGs for total IgE were identified in subjects with asthma (n = 1179 DEGs) than in all subjects (n = 631 DEGs). In subjects with asthma, DEGs were mapped to 11 gene modules. The top module for positive correlation with total IgE was mapped to myoepithelial and mucus secretory cells in lower airway epithelia and was regulated by IL-4, IL5, IL-13, and IL-33. Within this module, hub genes included CDH26, FETUB, NTRK2, CCBL1, CST1, and CST2. Furthermore, an enrichment analysis showed overrepresentation of genes in signaling pathways for synaptogenesis, IL-13, and ferroptosis, supporting interactions between interleukin- and acetylcholine-induced responses. CONCLUSIONS: Our findings for nasal epithelial gene expression support neuroimmune coregulation of total IgE in youth with asthma.

A novel piwi-interacting RNA associates with type 2-high asthma phenotypes.

BACKGROUND: Piwi-interacting RNAs (piRNAs), comprising the largest noncoding RNA group, regulate transcriptional processes. Whether piRNAs are associated with type 2 (T2)-high asthma is unknown. OBJECTIVE: We sought to investigate the association between piRNAs and T2-high asthma in childhood asthma. METHODS: We sequenced plasma samples from 462 subjects in the Childhood Asthma Management Program (CAMP) as the discovery cohort and 1165 subjects in the Genetics of Asthma in Costa Rica Study (GACRS) as a replication cohort. Sequencing reads were filtered first, and piRNA reads were annotated and normalized. Linear regression was used for the association analysis of piRNAs and peripheral blood eosinophil count, total serum IgE level, and long-term asthma exacerbation in children with asthma. Mediation analysis was performed to investigate the effect direction. We then ascertained if the circulating piRNAs were present in asthmatic airway epithelial cells in a Gene Expression Omnibus (GEO; www.ncbi.nlm.nih.gov/geo) public data set. RESULTS: Fifteen piRNAs were significantly associated with eosinophil count in CAMP (P ≤ .05), and 3 were successfully replicated in GACRS. Eleven piRNAs were associated with total IgE in CAMP, and one of these was replicated in GACRS. All 22 significant piRNAs were identified in epithelial cells in vitro, and 6 of these were differentially expressed between subjects with asthma and healthy controls. Fourteen piRNAs were associated with long-term asthma exacerbation, and effect of piRNAs on long-term asthma exacerbation are mediated through eosinophil count and serum IgE level. CONCLUSION: piRNAs are associated with peripheral blood eosinophils and total serum IgE in childhood asthma and may play important roles in T2-high asthma.

Health literacy and asthma: An update.

The US Department of Health and Human Services has defined health literacy (HL) as the degree to which individuals have the capacity to obtain, process, and understand basic health information and services needed to make appropriate health decisions. Structural and social determinants of health lead to low HL in approximately 36% of adults in the United States, where this condition is most prevalent in racial and ethnic minorities, economically disadvantaged communities, and immigrants with limited English proficiency. In turn, low HL can worsen asthma outcomes through direct effects (eg, nonadherence to or incorrect use of medications) and indirect effects (eg, an unhealthy diet leading to obesity, a risk factor for asthma morbidity). The purpose of this update is to examine evidence from studies on low HL and health and asthma outcomes published in the last 12 years, identify approaches to improve HL and reduce health disparities in asthma, and discuss future directions for research in this area under the conceptual framework of a socioecological model that illustrates the multifactorial and interconnected complexity of this public health issue at different levels.

Serum α-Klotho level, lung function, airflow obstruction and inflammatory markers in US adults.

Background: α-Klotho is a pleiotropic protein that may have anti-oxidative and anti-inflammatory properties in the lung, but its role in airflow obstruction or lung function is largely unknown. Methods: This was a cross-sectional study of 6046 adults aged 40-79 years in the US National Health and Nutrition Examination Survey (NHANES) 2007-2012. We used multivariable logistic or linear regression to examine the relation between serum α-Klotho level and airflow obstruction, defined as forced expiratory volume in 1 s (FEV1) <80% of predicted and FEV1/forced vital capacity (FVC) ratio <0.70; FEV1, FVC and FEV1/FVC as percentage of predicted; and inflammatory markers in blood (white blood cell count, eosinophils, neutrophils and C-reactive protein (CRP)). Results: α-Klotho levels in the second to fourth quartiles (Q2-Q4) were associated with significantly decreased odds of airflow obstruction (adjusted OR for Q2-Q4 versus lowest quartile (Q1) 0.54 (95% CI 0.35-0.81)) in never-smokers and ex-smokers with <10 pack-years of smoking, but not in current smokers or ex-smokers with ≥10 pack-years of smoking. In all participants, each unit increment in log10-transformed α-Klotho level was significantly associated with 5.0% higher FEV1 % pred and 3.7% higher FVC % pred. Higher α-Klotho was also associated with lower eosinophils, neutrophils and CRP in participants both with and without airflow obstruction. Conclusions: Higher serum α-Klotho is associated with lower inflammatory markers and higher lung function in adults with and without airflow obstruction, and with decreased odds of airflow obstruction in never-smokers and ex-smokers with <10 pack-years of smoking. Further studies are warranted to replicate our findings and evaluate underlying mechanisms.

Methyl-TWAS: A powerful method for in silico transcriptome-wide association studies (TWAS) using long-range DNA methylation.

In silico transcriptome-wide association studies (TWAS) are commonly used to test whether expression of specific genes is linked to a complex trait. However, genotype-based in silico TWAS such as PrediXcan, exhibit low prediction accuracy for a majority of genes because genotypic data lack tissue- and disease-specificity and are not affected by the environment. Because methylation is tissue-specific and, like gene expression, can be modified by environment or disease status, methylation should predict gene expression with more accuracy than SNPs. Therefore, we propose Methyl-TWAS, the first approach that utilizes long-range methylation markers to impute gene expression for in silico TWAS through penalized regression. Methyl-TWAS 1) predicts epigenetically regulated/associated expression (eGReX), which incorporates tissue-specific expression and both genetically- (GReX) and environmentally-regulated expression to identify differentially expressed genes (DEGs) that could not be identified by genotype-based methods; and 2) incorporates both cis- and trans- CpGs, including various regulatory regions to identify DEGs that would be missed using cis- methylation only. Methyl-TWAS outperforms PrediXcan and two other methods in imputing gene expression in the nasal epithelium, particularly for immunity-related genes and DEGs in atopic asthma. Methyl-TWAS identified 3,681 (85.2%) of the 4,316 DEGs identified in a previous TWAS of atopic asthma using measured expression, while PrediXcan could not identify any gene. Methyl-TWAS also outperforms PrediXcan for expression imputation as well as in silico TWAS in white blood cells. Methyl-TWAS is a valuable tool for in silico TWAS, leveraging a growing body of publicly available genome-wide DNA methylation data for a variety of human tissues.

Adverse Events during Adulthood, Child Maltreatment, and Asthma among British Adults in the UK Biobank.

Rationale: Intimate partner violence and child maltreatment have been separately associated with asthma in adults. No study has concurrently examined of adulthood adverse events (including, but not limited to, intimate partner violence) and child maltreatment on asthma in adults. Objectives: To concurrently examine of adulthood adverse events and child maltreatment on asthma in adults. Methods: This was a cross-sectional study of adulthood adverse events and child maltreatment on current asthma in 87,891 adults 40-69 years old who participated in the UK Biobank. Adulthood adverse events were assessed using questions adapted from a national crime survey. Child maltreatment was ascertained using the Childhood Trauma Screener questionnaire. Current asthma was defined as physician-diagnosed asthma and current wheeze and was further classified as noneosinophilic or eosinophilic according to eosinophil count (<300 vs. ⩾300 cells per microliter). Results: In a multivariable analysis, participants who reported two or more types of adulthood adverse events had 1.19-1.45 times significantly higher odds of asthma than those who did not, whereas participants who reported two or more types of child maltreatment had 1.25-1.59 significantly higher odds of asthma than those who reported no child maltreatment. After stratification by sex, similar results were obtained for child maltreatment in women and men, whereas adulthood adverse events were only significantly associated with asthma in women. Similar findings were observed in analyses that were restricted to never-smokers and former smokers with <10 pack-years of smoking and in analyses of noneosinophilic and eosinophilic asthma. Conclusions: In a cohort of British adults, child maltreatment was associated with current asthma in men and women, whereas adulthood adverse events were associated with current asthma in women only. This was independent of cigarette smoking or eosinophil count.

Epigenome-wide association analysis of infant bronchiolitis severity: a multicenter prospective cohort study.

Bronchiolitis is the most common lower respiratory infection in infants, yet its pathobiology remains unclear. Here we present blood DNA methylation data from 625 infants hospitalized with bronchiolitis in a 17-center prospective study, and associate them with disease severity. We investigate differentially methylated CpGs (DMCs) for disease severity. We characterize the DMCs based on their association with cell and tissues types, biological pathways, and gene expression. Lastly, we also examine the relationships of severity-related DMCs with respiratory and immune traits in independent cohorts. We identify 33 DMCs associated with severity. These DMCs are differentially methylated in blood immune cells. These DMCs are also significantly enriched in multiple tissues (e.g., lung) and cells (e.g., small airway epithelial cells), and biological pathways (e.g., interleukin-1-mediated signaling). Additionally, these DMCs are associated with respiratory and immune traits (e.g., asthma, lung function, IgE levels). Our study suggests the role of DNA methylation in bronchiolitis severity.

Distinction between rhinitis alone and rhinitis with asthma using interactomics.

The concept of "one-airway-one-disease", coined over 20 years ago, may be an over-simplification of the links between allergic diseases. Genomic studies suggest that rhinitis alone and rhinitis with asthma are operated by distinct pathways. In this MeDALL (Mechanisms of the Development of Allergy) study, we leveraged the information of the human interactome to distinguish the molecular mechanisms associated with two phenotypes of allergic rhinitis: rhinitis alone and rhinitis in multimorbidity with asthma. We observed significant differences in the topology of the interactomes and in the pathways associated to each phenotype. In rhinitis alone, identified pathways included cell cycle, cytokine signalling, developmental biology, immune system, metabolism of proteins and signal transduction. In rhinitis and asthma multimorbidity, most pathways were related to signal transduction. The remaining few were related to cytokine signalling, immune system or developmental biology. Toll-like receptors and IL-17-mediated signalling were identified in rhinitis alone, while IL-33 was identified in rhinitis in multimorbidity. On the other hand, few pathways were associated with both phenotypes, most being associated with signal transduction pathways including estrogen-stimulated signalling. The only immune system pathway was FceRI-mediated MAPK activation. In conclusion, our findings suggest that rhinitis alone and rhinitis and asthma multimorbidity should be considered as two distinct diseases.

Circulating MicroRNAs associated with Bronchodilator Response in Childhood Asthma.

Rationale: Bronchodilator response (BDR) is a measure of improvement in airway smooth muscle tone, inhibition of liquid accumulation and mucus section into the lumen in response to short-acting beta-2 agonists that varies among asthmatic patients. MicroRNAs (miRNAs) are well-known post-translational regulators. Identifying miRNAs associated with BDR could lead to a better understanding of the underlying complex pathophysiology. Objective: The purpose of this study is to identify circulating miRNAs associated with bronchodilator response in asthma and decipher possible mechanism of bronchodilator response variation. Methods: We used available small RNA sequencing on blood serum from 1,134 asthmatic children aged 6 to 14 years who participated in the Genetics of Asthma in Costa Rica Study (GACRS). We filtered the participants into high and low bronchodilator response (BDR) quartiles and used DeSeq2 to identify miRNAs with differential expression (DE) in high (N= 277) vs low (N= 278) BDR group. Replication was carried out in the Leukotriene modifier Or Corticosteroids or Corticosteroid-Salmeterol trial (LOCCS), an adult asthma cohort. The putative target genes of DE miRNAs were identified, and pathway enrichment analysis was performed. Results: We identified 10 down-regulated miRNAs having odds ratios (OR) between 0.37 and 0.76 for a doubling of miRNA counts and one up-regulated miRNA (OR=2.26) between high and low BDR group. These were assessed for replication in the LOCCS cohort, where two miRNAs (miR-200b-3p and miR-1246) were associated. Further, functional annotation of 11 DE miRNAs were performed as well as of two replicated miRs. Target genes of these miRs were enriched in regulation of cholesterol biosynthesis by SREBPs, ESR-mediated signaling, G1/S transition, RHO GTPase cycle, and signaling by TGFB family pathways. Conclusion: MiRNAs miR-1246 and miR-200b-3p are associated with both childhood and adult asthma BDR. Our findings add to the growing body of evidence that miRNAs play a significant role in the difference of asthma treatment response among patients as it points to genomic regulatory machinery underlying difference in bronchodilator response among patients. Trial registration: LOCCS cohort [ClinicalTrials.gov number: NCT00156819], GACRS cohort [ClinicalTrials.gov number: NCT00021840].

Asthma and Chronic Obstructive Pulmonary Disease.

In the United States, asthma and chronic obstructive pulmonary disease (COPD) disproportionately affect African Americans, Puerto Ricans, and other minority groups. Compared with non-Hispanic whites, minorities have been marginalized and more frequently exposed to environmental risk factors such as tobacco smoke and outdoor and indoor pollutants. Such divergent environmental exposures, alone or interacting with heredity, lead to disparities in the prevalence, morbidity, and mortality of asthma and COPD, which are worsened by lack of access to health care. In this article, we review the burden and risk factors for racial or ethnic disparities in asthma and COPD and discuss future directions in this field.

Association of World Trade Center (WTC) Occupational Exposure Intensity with Chronic Obstructive Pulmonary Disease (COPD) and Asthma COPD Overlap (ACO).

INTRODUCTION: Reported associations between World Trade Center (WTC) occupational exposure and chronic obstructive pulmonary disease (COPD) or asthma COPD overlap (ACO) have been inconsistent. Using spirometric case definitions, we examined that association in the largest WTC occupational surveillance cohort. METHODS: We examined the relation between early arrival at the 2001 WTC disaster site (when dust and fumes exposures were most intense) and COPD and ACO in workers with at least one good quality spirometry with bronchodilator response testing between 2002 and 2019, and no physician-diagnosed COPD before 9/11/2001. COPD was defined spirometrically as fixed airflow obstruction and ACO as airflow obstruction plus an increase of ≥ 400 ml in FEV1 after bronchodilator administration. We used a nested 1:4 case-control design matching on age, sex and height using incidence density sampling. RESULTS: Of the 17,928 study participants, most were male (85.3%) and overweight or obese (84.9%). Further, 504 (2.8%) and 244 (1.4%) study participants met the COPD and ACO spirometric case definitions, respectively. In multivariable analyses adjusted for smoking, occupation, cohort entry period, high peripheral blood eosinophil count and other covariates, early arrival at the WTC site was associated with both COPD (adjusted odds ratio [ORadj] = 1.34, 95% confidence interval [CI] 1.01-1.78) and ACO (ORadj = 1.55, 95%CI 1.04-2.32). CONCLUSION: In this cohort of WTC workers, WTC exposure intensity was associated with spirometrically defined COPD and ACO. Our findings suggest that early arrival to the WTC site is a risk factor for the development of COPD or of fixed airway obstruction in workers with pre-existing asthma.

Systems Genomics Reveals microRNA Regulation of ICS Response in Childhood Asthma.

BACKGROUND: Asthmatic patients' responses to inhaled corticosteroids (ICS) are variable and difficult to quantify. We have previously defined a Cross-sectional Asthma STEroid Response (CASTER) measure of ICS response. MicroRNAs (miRNAs) have shown strong effects on asthma and inflammatory processes. OBJECTIVE: The purpose of this study was to identify key associations between circulating miRNAs and ICS response in childhood asthma. METHODS: Small RNA sequencing in peripheral blood serum from 580 children with asthma on ICS treatment from The Genetics of Asthma in Costa Rica Study (GACRS) was used to identify miRNAs associated with ICS response using generalized linear models. Replication was conducted in children on ICS from the Childhood Asthma Management Program (CAMP) cohort. The association between replicated miRNAs and the transcriptome of lymphoblastoid cell lines in response to a glucocorticoid was assessed. RESULTS: The association study on the GACRS cohort identified 36 miRNAs associated with ICS response at 10% false discovery rate (FDR), three of which (miR-28-5p, miR-339-3p, and miR-432-5p) were in the same direction of effect and significant in the CAMP replication cohort. In addition, in vitro steroid response lymphoblastoid gene expression analysis revealed 22 dexamethasone responsive genes were significantly associated with three replicated miRNAs. Furthermore, Weighted Gene Co-expression Network Analysis (WGCNA) revealed a significant association between miR-339-3p and two modules (black and magenta) of genes associated with immune response and inflammation pathways. CONCLUSION: This study highlighted significant association between circulating miRNAs miR-28-5p, miR-339-3p, and miR-432-5p and ICS response. miR-339-3p may be involved in immune dysregulation, which leads to a poor response to ICS treatment.

Cis- and trans-eQTM analysis reveals novel epigenetic and transcriptomic immune markers of atopic asthma in airway epithelium.

BACKGROUND: Expression quantitative trait methylation (eQTM) analyses uncover associations between DNA methylation markers and gene expression. Most eQTM analyses of complex diseases have focused on cis-eQTM pairs (within 1 megabase). OBJECTIVES: This study sought to identify cis- and trans-methylation markers associated with gene expression in airway epithelium from youth with and without atopic asthma. METHODS: In this study, the investigators conducted both cis- and trans-eQTM analyses in nasal (airway) epithelial samples from 158 Puerto Rican youth with atopic asthma and 100 control subjects without atopy or asthma. The investigators then attempted to replicate their findings in nasal epithelial samples from 2 studies of children, while also examining whether their results in nasal epithelium overlap with those from an eQTM analysis in white blood cells from the Puerto Rican subjects. RESULTS: This study identified 9,108 cis-eQTM pairs and 2,131,500 trans-eQTM pairs. Trans-associations were significantly enriched for transcription factor and microRNA target genes. Furthermore, significant cytosine-phosphate-guanine sites (CpGs) were differentially methylated in atopic asthma and significant genes were enriched for genes differentially expressed in atopic asthma. In this study, 50.7% to 62.6% of cis- and trans-eQTM pairs identified in Puerto Rican youth were replicated in 2 smaller cohorts at false discovery rate-adjusted P < .1. Replicated genes in the trans-eQTM analysis included biologically plausible asthma-susceptibility genes (eg, HDC, NLRP3, ITGAE, CDH26, and CST1) and are enriched in immune pathways. CONCLUSIONS: Studying both cis- and trans-epigenetic regulation of airway epithelial gene expression can identify potential causal and regulatory pathways or networks for childhood asthma. Trans-eQTM CpGs may regulate gene expression in airway epithelium through effects on transcription factor and microRNA target genes.